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1-^ I Abstract: The asymmetric exclusion process on a ring in one-dimension is considered with 

^1 ■ a single defect particle. The steady state has previously been solved by a matrix product 

method. Here we use the Bethe ansatz to solve exactly for the long time limit behaviour 

of the generating function of the distance travelled by the defect particle. This allows us to 

recover steady state properties known from the matrix approach such as the velocity, and 

^ I obtain new results such as the diffusion constant of the defect particle. In the case where the 

^r^ ■ defect particle is a second class particle we determine the large deviation function and show 



that in a certain range the distribution of the distance travelled about the mean is Gaussian. 

Moreover the variance (diffusion constant) grows as L^/^ where L is the system size. This 
G^ i behaviour can be related to the superdiffusive spreading of excess mass fluctuations on an 

^ ' infinite system. In the case where the defect particle produces a shock, our expressions 

c^ ■ for the velocity and the diffusion constant coincide with those calculated previously for an 

infinite system by Ferrari and Pontes. 
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1. Introduction and Model Definition 

The asymmetric simple exclusion process (ASEP) [|l| is a simple example of a driven 
lattice gas and as such is a system far from thermal equilibrium. The model comprises 
particles hopping in a preferred direction along a one dimensional lattice with hard core 
exclusion imposed. The model's broad interest lies in its connections to growth processes, 
the problem of directed polymers in a random medium and Burgers equation P, §. It is 
also a natural starting point for many traffic flow models . 

Multi-species variants of the ASEP have been considered [^ |^, ^, |^, 0. In particular the 
idea of a second class particle has proven useful [^]. The second class particle hops forward as 
usual when the neighbouring site is empty but is overtaken by the other particles. Therefore 
it moves forward in an environment of low density of particles and backwards in a high 
density environment. In this way a second class particle can be used to locate shocks which 
are sudden changes in density over a microscopic region [jll], |T^, |T3|, |T^. A generalisation 



of the second class particle idea to that of a defect particle |T^, |16] has been shown to 
exhibit phase transitions and, in particular, phase coexistence. Interpreted in the context of 
traffic problems, the phase transition corresponds to the appearance of a traffic jam whereas 
coexistence between phases of different densities corresponds to the coexistence between a 
freely flowing and a jammed region in traffic. The model has also been interpreted in the 
context of a two way road 



Analytical results for the ASEP, such as the diffusion constant |18[, and for the second 



class particle problem such as the steady state current |P, |T5|, |T6l have been obtained via a 



matrix product technique [0, Q. Recently it was shown that for the ASEP, all moments 
and the large deviation function of the time integrated current can be computed via a Bethe 
ansatz technique |^|. Here we show that Bethe ansatz can also be used in the case of the 



defect particle and this allows us to generalize previous ||2^, |TB|, |TB| results obtained by the 
matrix approach. 



Let us now define the model we consider |T5|, |T^. The model comprises a single defect 
particle (indicated by 2) and M — 1 first class {i.e. normal) particles (indicated by 1) on a 
ring of size L sites. The hopping rates of the particles are as follows 

10^01 with rate 1 
2 0^02 with rate a 
12^21 with rate /? . (1) 

By this it is implied, for example, that in an infinitessimal time interval dt a first class 
particle hops to the neighbouring site to the right with probability dt if that neighbouring 
site is empty. We restrict ourselves to a, /3 > 0. 

1 



2. Main Results 



Before discussing the technical details of the Bethe ansatz solution we summarise in this 
section our main results. Let us denote by yt the distance travelled (total number of hops 
forward minus total number of hops backward) by the defect particle. In the steady state, 
yt is a random variable. Its first and second moments give the velocity v and the diffusion 
constant A of the defect particle 
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All cumulants of yt can be computed from the knowledge of the generating function {e"'^^) 
via 
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Here, by employing a Bethe ansatz technique we calculate exactly the large t behaviour 
of this generating function namely 



^'^ J->oo t 
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for arbitrary L and M. Exact expressions of v and A follow easily from the knowledge of 

A(7) 
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In the thermodynamic limit {L, M large) with fixed density p where 

p = M/L 



(7) 



our exact expressions for A (7) allow us to show that the velocity v and diffusion constant 
A of the defect particle have the following asymptotic forms in different regions of the 
parameter space oi a, P, p 



For (3 > p> 1 - a 
For P < p and 1 — a < p 
For P > p and 1 — a > p 
For P < p <l - a 



V = 1 — 2p and A 

V = 1 — P — p and A 

V = a — p and A 

V = a — P and A 
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These results lead to a phase diagram which is essentially the same as the one given 
in 



[p!q] with the same expressions for the velocity. For /?>p>l — aitis known p, |T6 
that the density profile, as seen from the defect particle, has a power law decay towards its 
asymptotic value. In the phases where P,l — a > p and (3,1 — a < p the density profile 
decays exponentially towards the asymptotic values. In the final phase there is coexistence 
between a region of low density /? in front of the defect particle and high density 1 — a behind 
the defect particle. Therefore there is a shock separating the two regions at a distance xL 



in front of the defect particle where x is given by p = /5x + (1 — a)x [^, |T6[. The novelty of 
the present work is that, as we can calculate A(7), all the higher cumulants of the distance 
travelled (including the diffusion constant) can be calculated exactly for the different phases. 

One should notice from (^) that in the whole phase (3 > p and l—a<p (which includes 
the case of a second class particle (3 = a = 1) the diffusion constant of the defect particle 
increases with L. This is in contrast with the diffusion constant of the first class particles 
|T8| , 21| which (in the absence of a defect) decreases as L"^'"^. This difference between first 



and second class particles is not a surprise since, for an infinite system, one expects the 
fluctuations of position to be super diffusive for second class particles and subdiffusive for 
first class particles [^ ^. 

Furthermore, in the whole phase 1 — a < p < j3 we will show that A (7) defined by (^ is 
given in the large L limit by 



A(7)-7(l-2p):^^ 



on the scale where 'y ^ 0{L 
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The fact that in a certain range 7, the expression of A(7) is quadratic implies that the 
distribution of the variable yt/t — (1 — 2p) is Gaussian over a certain range. It is easy to 
check that the range which corresponds to 7 ~ 0(L^^/^) is yt/t — (1 — 2p) ~ 0(1/L), so 
that over that range, the distribution of this difference should be Gaussian. 

In the following sections we give the derivation of (PHT^). 



3. Generating function for fluctuations in distance travelled 

To calculate A (7), we follow and extend the technique of references PTl ^^. First consider 



Pt{C,y) which is the probability of the system being in configuration C and of the defect 
particle having been displaced a distance y (negative if the particle has been displaced 
backwards). The master equation is 

dPt{C,y) 
dt 

Y: [Mo{C , C')P,{C', y)+M,{C, C')P,{C\ y-l)+ M^,{C , C')P,{C', y + 1)] (13) 
c 

where J^o{C,C'), M.i{C,C'), M.-i{C ,C') are rates for transitions from C to C that do 
not move the defect particle, move the defect particle forward, move the defect particle 
backward respectively. The total rate out of configuration C is given by — A^o(C,C) = 
EcYc[-^i('^'.'^) + Mo{C',C) + M-i{C',C)]- The generating function 

+00 
F,{C)= Y^ exp(72/)Pt(C,2/) (14) 

y=-oo 

obeys 

^^j^ = Y [Mo{C , C')F,{C') + e^MiiC , C')F,{C') + e'^M-iiC , C')F,{C')] . (15) 

Now 

(e^^') = ^F,(C) (16) 

c 

so we expect as in (|^) 

where A(7) is the largest eigenvalue of the matrix A4{j) = A^o + e'^-Mi + e~'^M.^i. For 
7 = we know that A(7) = because A^(0) is a stochastic matrix. Now as in 1^1], ^, by the 



Perron- Frobenius theorem we know that the largest eigenvalue of M. (7) is non-degenerate 
therefore as 7 increases from zero there can be no crossing of the largest eigenvalue. Thus 
A(7) is the eigenvalue that tends to zero as 7 tends to zero. 

One should note in the present problem that for 7 = the eigenvector with eigenvalue 
zero (the steady state) is non-trivial and has previously been constructed by using a matrix 
product [^, |16|. In the following sections we show how the Bethe ansatz can recover some 
properties of this steady state in the limit 7 — > 0. 



4. Bethe Ansatz 

Let a configuration of the particles on tlie ring be specified by the co-ordinates {xi, X2, ■ ■ ■ , xm} 
where xi is the position of the defect particle and X2 < x^ . . . < xm are the positions of the 
normal particles. By convention, one can always choose {xi . . .xm} such that 1 < xi < L 
and Xi < X2 < . . . < Xm < Xi + L. With the dynamics (|T]) the equation for an eigen- 
function ip{xi, . . . xm) of ^^(7) with eigenvalue A (7) takes the form: for Xi < Xj+i — 1 and 
Xm < L + Xi — 1 (i.e. when all particles are more than one lattice spacing apart) 

M 

\i){xi, . . . Xm) = -{M-l+a)ij{xi, . . . XM)+e'(xi'{xi-l, . . . XM)+^i'{xi, ..., x^-l, . . . xm) 

(18) 
When Xi = Xj+i — 1 or xm = Xi + L — 1 (i.e. for configurations where two consecutive sites 



are occupied), the equation (18) is in principle modified. Insisting that it remains valid even 
in these particular cases, requires that the function ipi^xi, . . . xm) take values in unphysical 
regions (xj+i = Xi or xm = Xi + L) which satisfy the following conditions arising from the 
interaction of particles: 

(1 - /5)z^(xi, ...xi + L-l)- e^aip{xi - 1, . . . xi + L - 1) = (19) 

^{...,Xi,Xi + l,...)-i){...,Xi,Xi,...) = Q for l<i<M (20) 

ail){xi,xi + 1 . . .) - ^(xi, xi . . .) + e~^l3il){xi + 1, X3, . . . , Xi + L) = (21) 

The Bethe ansatz consists of writing the eigenfunction V^(xi, . . . xm) as 

ij{x,,...XM)=e^'''a'''Y.Mi)-QiM){zQ^i)T---{zQ(^M)r" . (22) 

Q 

where the sum is over all permutations Q of 1 . . . M. The amplitudes Aq(i)...q(^m) and the 
wave numbers Zi. . . zm are a priori arbitrary complex numbers. This ansatz inserted into 



(p!8|) gives for the eigenvalue 

M 1 

A(7) = -(M-1 + «) + $:-. (23) 

fc=i ^^ 

for any choice of the amplitudes ^q(i)...q(a/) and of the wave numbers Zi . . . Zm- For (|^) to 
fulfil conditions (p!9|-plD, the amplitudes and the wave numbers have to satisfy 

z^ [(1-/3)^,-1] 
z^ [(l-/5)z,-l] 



•^i---« ~ ^ ' L {(-[ o\ 11 "^i-j \ I 



Zj-1 



J\...ji... — I ) J\....ij... l'^"3) 

Zi — l 

Aji... = (— )- — - {azj — l)Aij... + aPzizfAi...j + aPzjZ^Aj...i . (26) 
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Using ([2^ . |25|) allows ( ^61) to be written as 

l + a/3 
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Using ( [27| ) twice in succession yields, after some algebra, the following condition on the 
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Af-l 



M 



a 



PUi^-Zk] 



k=l 



1 — Zi 



{azj - 1) {bzj -1) {zj-1) 



M-l 



M 
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^J k=l 

The wavefunction '?/'(xi, . . . xm) corresponding to the largest eigenvalue A(7) is invariant 
under translation, therefore 

M 

e''al[zk = l. (30) 



k=l 



One can rewrite the Bethe equations (p9|) in terms of two constants C, E as follows 
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(31) 
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Under this form the Bethe equations (^) are much easier to solve. One first finds the Zi 
solutions of (|3^) (which depend on the unknown constants C and E). Then by inserting 
these solutions into (|30D and (0), the constants C and E are determined. In Appendix A, 
we show ( |A14| , |Al5|jA19D that, in so doing, the eigenvalue A (7) can be written as 



A(7) = 


00 r^n 
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a _ 


^. ^l«^'l' 
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,_ -z'^[l + {z-l)E] 



where 

^^^^" {bz-l){az-l){z-lY'-^ 
and the constant E is fixed by imposing 

^ C" r / / 1 d^ 1 
^ n /i /i 27rz z 



-y[Q(^) 



(33) 
(34) 

(35) 
(36) 



As is shown in the appendix A the contours of integration in (|33| , |3^ , |36D are small contours 
which surround 1 and 1/a but do not surround 1/6 ; the particular cases where a = 1, b = 1 
or a = b can be obtained easily as limiting cases since all the integrals which appear in the 
right hand side of (p3|,pl|,p6D are rational functions of a and /3. 



Equations (^,34,35,^) determine the exact expression of A(7) for arbitrary L, M, a and 
f3. The difference between these equations and the corresponding equations of |^| is that 
in the present case we have an additional unknown constant E. This feature emerges from 
the structure of the Bethe equations (P^). 



5. Exact Expressions for the Velocity and Diffusion Constant 



In principle, one can use (^) to expand E in powers of C. Replacing E by its expansion 
in powers of C in (^,0) gives the expansions of A and 7 in powers of C. Then one can 
eliminate C between the two expansions and this gives A in powers of 7. This is what is 
done in this section to obtain exact expressions of the velocity and of the diffusion constant. 



For example to obtain the velocity f , one can note from ( P^ that C vanishes linearly 
with 7 so that from (|36|) , the limiting value E{0) of -E at 7 = is 



where Xl^m is defined by 
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(38) 



27ii {z ~ ly^ {bz - l){az - 1) ' 

Then from (|33| , |3^ , one finds that for 7 small A (7) = f 7 + 0(7^) with the velocity v given 
by 

Xl,mXl-2,M-2 — Xl M-iXl-2,M-'^ /on\ 

V = — ■ ■ ■ ■ (39) 



'L,M 



where 



Expression 



to obtain 



Zl,m — XlmXl^i^M-2 — Xl^m-iXl-i,m-i ■ 
may be simplified by using 

Xl,M = Xl^im + -^L-l,A/-l 

Zl-i,m — Zl-i.m-i 

v = ^ ^ — • 



(40) 



(41) 



(42) 



In a similar fashion one obtains from the second derivatives of ( ^PH^Gl) at 7 = 0, after a 
good deal of straightforward but tedious algebra, 
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(43) 



where W2l,2M is defined by 

W2L,2M = 



+ 



dz 



.2L 



[i+Emz-i)f 



2m {z - If^ {hz - lY{az - l)'^ 
Expression (^) can be simplified by using (|39|^ID to obtain 
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(45) 



Alternatively one could obtain (^51) directly from ( p3l - pq ) by using, for example, Mathemat- 
ica. The integrals in the above expressions can be evaluated by residues. In this way one 
can show that the integral expressions for Z^^m and v given by (E3,En,H^) are equivalent to 



those derived in [|T5|, |T6| within the matrix product formulation. However exact evaluation 
of the integrals involved in the diffusion constant (^5|) results, in general, in cumbersome 
expressions. 



Remark: For the case of a second class particle {a = (3 = 1) simplification of 
possible and one recovers the expression first presented in p2|: 
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(2L-3)! 



(2M-1)!(2L-2M+1)! 



(M-1)!(L-M)! 



X [{L - 5)(M - 1)(L - M) + (L - 1)(2L - 1)] . 



(46) 



However the derivation of (|B) from ( ^3] ) is tedious and not illuminating therefore we do not 
present it. 



6. Asymptotics and Phase Diagram 

In order to obtain tlie ptiase diagram it suffices to determine the asymptotic form of Z^^m 
which in turn determines the asymptotic form of the velocity (^). We carry this out in 
detail in appendix B, but would like to outline here how the different phases arise. We 
restrict ourselves to /? < 1 (6 > 0) and a < 1. 

Consider the quantity Xl^m given by (^). In the limit of L^M large (with fixed p = 
L/M) there are three possible dominant contributions to the integral, all lying on the real 
axis: a saddle point at Zc = 1/(1 — p); a pole aX z = 1/6 and a pole at z = l/a. The possible 
dominant and subdominant contributions to X^^m gives rise to four phases as follows. 

1. Ill/ a < Zc < l/h the contours of the two integrals in ( pS]) may be merged and deformed 
to pass through the saddle point. Therefore only the saddle point contributes. 

2. If 1/a < Zf. and 1/6 < z^. the contours of the two integrals may be merged and deformed 
to pass through the saddle point. However the contour must make a clockwise detour 
around the pole at 2; = 1/6. Therefore the pole at 2; = 1/6 is the dominant contribution 
and the saddle point is the subdominant contribution. 

3. If 1/a > Zc and 1/6 > z^ the two integrals give separate contributions. The pole 
sX z = 1/a gives the dominant contribution and the integral around z = 1 may be 
deformed to pass through the saddle point and gives the subdominant contribution. 

4. If 1/a > Zc > 1/6 the two integrals give separate contributions: the pole ai z = 1/a 
and the clockwise integral around the pole at 2; = 1/6. 

In appendix B the dominant contributions to the desired integrals are evaluated and 
expressions (§]-|TT|) are established. At this point we can already see the interesting feature 
that since Z in (^Dp is a difference of products of X, the subdominant contribution as well as 
the dominant contribution to the integral X must be evaluated to obtain Z and the velocity. 
Also note that in phase where 1/a < z^ < 1/6 , power law decays in correlation functions, 
for example the density profile, will arise from the saddle point being dominant. In the other 
phases a dominant pole will give rise to exponential decays. Of particular interest is the 
phase where 1/a > 2;c > 1/6 and the two poles compete. As described in the introduction 
this is the phase where a shock exists. 



7. Scaling of the Large Deviation Function for 1 — a < p < (3 



In this phase, because the integrals are dominated by the saddle point, the analysis of 
the asymptotics of the exact expression of A (7) given by (|33|-|36D is rather different from the 
other cases. In these expressions a contour integral actually implies two integrals around 
z = 1 and z = 1/a but in this phase, for large L, one expects all the integrals to be 
dominated by their saddle point Zc = 1/(1 — p)- So to lighten the notation we write a single 
integral. Let us replace the variables z and E in (B3-p^ by y and F 



z = Zc + y 



(47) 



and 



E 



1 1 F 

+ 



Zc — 1 Zc — 1 L 

Clearly, the values of y which contribute to the integrals in (|3^j3^j3^) are of order y 
0{L~^^'^) so that one can rewrite (p3| , p4| , p6D as 



(48) 
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where Sp is given by 



s, = E 



C" r dy 
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n J 27Ti 



R{y) 



F ^ Fy 



with 



R{y) 



1 



Zc - 1 L L{zc-l)J 
{zc + yf 



(52) 



(53) 



;i - hzc - hy){azc + ay-l){zc + y- 1)*^-^ ' 
Under the assumption (which we will check later) that F is of order 1 (in the large L limit), 
if we define 

giy) = iog(^c + y)- p iog(zc - 1 + y) (54) 

we can evaluate the leading orders of the integrals ([5^ ) . 



For p odd, the leading large L behaviour is given by 



Sp — 



:-)^i-P D 



1 + 2 



'l-K p L 



1+f 



.^"(0), 



p\ 



(55) 
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where 



D = Ce^3(°)- 



1 - Zr. 



{bzc - l){azc - 1) 
For p even the leading order in the range where D is of order 1 is 



(56) 
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27rL^ \9"{0)J 
x|f/(0)2(p-1)!! + 



:i-pf 



a 



p p + b — 1 a — 1 + p 



g"{0){p + l)\\ 



1-P9'"i0) 
p 6 



(p + 3)!! 



D^fl-pY 1 



P+3 



P+2 



ip+l)\ 



2 V P / L^ VV(0) 



where we define {p — 1)!! = 1 for p = 0. F is fixed through (|5TD by 

5i ~ (Ze - 1)^0 



(57) 



and from (^ , pOD we obtain 
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4i^c-l 



-^1 



A(7)-7 



2- Zr 



zi{z,- ly 



-\S^{2-'iz,)-S2Z,{l-z,)]. 
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(59) 
(60) 



Using ( [55| ) and (|57D, ( ^8] ) gives 



(1-p) 



a 



1 + p 6 

p 0— 1+p a— 1+p 



+ 



D 



p(l - p)25/2 



which is consistent with our earher asumption that F is of order 1. From 



-D 



7 



A(7)-7(l-2p) 



27rp(l - p)L3/2 



(61) 

we find 

(62) 



27 
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2-p 1 
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p- 13 a -I + p 



^7rp(l-p) 



+ ^2^1/2 ^^J 1^ (63) 



as announced in (|T2D. 



11 



8. Discussion 

In this work we have shown how the Bethe ansatz can be used to calculate exactly via 
(p3|-|36D the large deviation function of the displacement of the defect particle. By analyzing 
the asymptotics of (|33| -p^) we could recover the velocity and phase diagram |l^, |T^ of the 
asymmetric exclusion model with a moving defect. The approach also allows new results 
(such as all the cumulants of the displacement of the defect) to be obtained, in particular 
the diffusion constant of the defect particle in the various phases. This adds to the body 
of knowledge concerning diffusion constants within the asymmetric exclusion process. For 
example the exact expression for the diffusion constant of a first-class particle calculated in 
||T8| allowed the determination of a universal amplitude for the centre of mass fluctuation 
for growth processes described by the one dimensional KPZ equation [§]. The diffusion 
constant of a second class particle is of interest since it is closely related to the motion of 



shocks. In |T^ the diffusion constant for a second class particle starting at the origin of an 
infinite lattice with a shock initial condition was calculated. Our results ( |Tll) for the phase 
exhibiting a shock {(3 < p < 1 — a) exactly agrees with that of \TA\. This is of interest 



since it shows that a single defect can provoke a shock in the ring geometry with the same 
behaviour as a shock on the infinite line pOf. If the fluctuations of the shock on a ring with 



a defect are identical to the fluctuations of the shock on an infinite line, this means that 



our results (|33| - |36D should give the whole large deviation function of a shock position on an 
infinite line. The behaviour of shocks and shock fluctuations is also connected to the phase 
diagram and density profile for systems with open boundary conditions 



When 7 — > the wave function ( ^4 ) reduces to the steady state probabilities for each 
configuration. Therefore, in principle, the steady state of the system, previously constructed 
using a matrix product p, |1^, |1^ , can be extracted from the present Bethe ansatz (p^) by 
taking the limit 7 — > 0. This shows that a co-ordinate Bethe ansatz is capable of describing 
non-trivial steady states of stochastic systems. In particular it would be interesting to 
understand further how this works and to see if the approach might be generalisable to 
larger numbers of species. The matrix product steady state of the present model is very 
closely related to that of the ASEP with open boundary conditions |1^, ^j. It would be 
of great interest to determine whether some implementation of the Bethe ansatz, perhaps 
related to that of the present work, could be used to recover the steady state with open 
boundary conditions. A major difficulty in doing so is that the particle number is not 
conserved with open boundaries. 

Returning to the case of a second class particle it is of interest to review how its dynamics 
are related to the spreading of excess mass. The central idea, termed coupling |l|], is well- 
known in the mathematical community but less so within physics. Consider two systems 

12 



containing only first class particles, identical except that one system has M particles and 
the other M — 1 particles. The two systems start from initial conditions differing only by 
the position of the extra particle in the system with M particles. In order to implement the 
dynamics one can consider at each time step randomly choosing a pair of sites i,i + 1 to 
update; then if there is particle at site i and a hole at site i + 1 the particle is moved forward. 
In the dynamics let us choose the same pairs of sites in the two systems at each update (one 
can think of using the same random numbers in a Monte Carlo program). Now it is easy 
to convince oneself that after any length of time the configurations of the two systems will 
differ only by the position of the extra particle (note that if we label the particles, the label 
of the extra particle will change under the dynamics). Further, one can convince oneself 
that the position of the extra particle has precisely the dynamics of a second class particle 
in the ASEP. Conversely, the system comprising M — 1 first class particles and one second 
class particle that we have studied describes the motion of an extra particle added to a 
system of M — 1 particles. Therefore the diffusion constant of the second class particle we 
have calculated here serves to describe the spreading of excess mass in the ASEP. 



Approximate calculations such as mode coupling |28, 23] have led to the following un- 
derstanding of the motion of excess mass fluctuations: the drift speed is 1 — 2p as in (^ for 
a = P = 1 and the spreading of density fluctutations around the drift grows as t^/^ on an 
infinite system i.e. it is superdiffusive. This superdiffusive motion can be recovered from 
the L^/^ finite system size dependence (|) of the diffusion constant of a second class particle 



2J] if we assume that a scaling form holds and the variance of the distance travelled by the 



second class particle can be written as 

(z/5 - {ytf - tL'"f{tlV) (64) 

where z is the dynamic exponent and f{x) is a scaling function tending to a constant as 
X ^ cxD. Now the dynamic exponent for the ASEP is known by the Bethe ansatz to be 
z = 3/2 [^, ^ and we expect the same exponent in the present model. In the limit L ^ oo 



for large but fixed t, the variance (y^) — {yt)"^ should not depend on system size therefore 
the scaling function must obey f{x) ~ x^^^^ as a; — > 0. We then find in this infinite system 
limit that (y^) — (yt)"^ ~ t^^'^ so that the typical spread of density fluctuations grows as t^/^. 
The spreading of mass fluctuations is also related to the scaling length ^ ~ t^/^ of the KPZ 
equation in one dimension (see |^ for detailed discussion). 

Finally let us mention that one can easily extend the calculation of the present paper to 
calculate the joint distribution of the distance yt covered by the defect particle and of the 
total distance Yt covered by all the first class particles. One can show that 



In 



/Q"/yt+SYt\ 



X{^,6)=\im^— '- (65) 

t-»oo t 
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is still given by (|23|) for arbitrary L and M where the Zi and the constants C and E satisfy 

M 



,7+(M-l)<5 



a Jl Zfc = 1 
fc=i 



C 



M+l ^ n JL 



M 



a 



/3e^^ n(^fc-i) 



(66) 



(67) 



fc=i 



instead of (|30| . |3T|) with (^) unchanged. 
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Appendix A: Analysis of the Bethe ansatz equations 



The solution for {zi} of (^0|- |32|) which gives A ^ as 7 ^ is of the form zi -^ 1/a and 
Zk —>■ 1 for 2 < k < M. We analyse separately the cases a 7^ 1 and a = 1 and for simplicity 
assume /? 7^ 1 and j3 ^ a although there is no problem in extending the analysis to include 
these cases. 



Case a ^ 1 
Consider the root zi and the M — 1 roots z/c for 2 < A; < M of 



which we rewrite as 



{z - lf'-\hz - l){az - 1) + z^C[l + {z- l)E] = 

such that Zi — ^ 1/a and 2;^ ^ 1 as C — > 0. Define 

^ ~z^[l + {z-l)E] 
^' {hz-l){az-l) ' 

then if z^ is a root of ( |A1| ) such that z^ — !► 1 as C — * one has for small C 

zu = l+ [Ce^'^'^Rizk) 

We wish to calculate expressions (^,^,0) of the form Y.k h{zk) (e.g. Eq. 
1/z). If h{z) is analytic near z = 1, one has by the residue theorem 



h{zk) 



dz 1 

— h{z)- 
1 27ri 



Ce^^^^R{zk) 



A/-1 r>/ 



R'{z)/[{M-l)R{z)] 



z-l-[Ce^^^^R{zk)]' 
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(Al) 



(A2) 



(A3) 
where h = 

(A4) 



where the contour is a circle centred on 1 and of radius e with \C\ ^-i ^ e ^ 1. Expanding 

1 
in powers of C '^^-^ and after an integration by parts, this gives 



"(-)-'W^£^[--1'^il^"'w^ 



p=iP 
Then summing over the roots 2 < k < M leads to 

E Kz,) = {M- 1)M1) + Y.-cH,r-h 

k=2 n=\ ^ "^1 ^'^^ 

Similarly, if h[z) is analytic near z = 1/a and we define 

-z^[l + {z-l)E] 



[Riz)r 



(A5) 



(^-1) 



n(M-l) 



(A6) 



Siz) 



then 



h{zi) 



a{bz - l){z - 1)^^-^ 



dz ^ ^ ^ 1 - CS'(z) 
h[z) 



(A7) 



(A8) 



h/a 2m ^ z-\- CS{z) 
where the contour is a circle centred on \ja and of radius e with |C| ^ e <^ 1. After 
expanding in powers of C and an integration by parts this gives 



(A9) 







Mzi) = h — 


CO ^n 


k 

a 


dz 
27ri 


h' 






Therefore, 


if one defines Qi^z) by 
















Q{z)- 


-z^[l + 


[z - l)E] 




T ~ 


(^ 


R{z) 


S{z) 




{hz-l){az- 


i)(^-i: 


M- 


~ _ i 

a 


one finds by combinin 


g (A|^A|) 















M 



Y^h{zk) = {M-l)h{l) + h 



k=l 



a 



CO /~in 



n=l 



n 



dz 
2Td 



h\z)[Q{z)\ 



(AlO) 



(All) 



Writing (H),® as 



one finds that 



M 



A(7) = -(M-l + «) + ^- 

fc=i ^fc 



M 



7 = - log a - XI log Zk 
fc=i 



oo f^n 

A(7) = -E- 

oo /-in 



+ 



2m z-^ 



n=\ 



n 



dz 1 
27ri z 



\Q{^)r 



(A12) 
(A13) 

(A14) 
(A15) 
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where Q{z) is given by (|A10|) . Then, with the use of (|A3|) , replacing (|3T|) by 



M 



1 



k=2 «/^ 



and using the fact ( [AB| , [A9D that 



M 



= \n{aP)+\n{l-z^) + J^-^J2^nR{zk) 

1 



k=2 



oo^ C*" 1 



,,^M&lu,i. ' I^Wl" 



f^ n 27ri [M - 1 Ii [z - 1)^^-^) Ji. z-l 
one finds that (|3T|) is satisfied if 

dz 1 



z-^Y 

a' 



o=i: 



n=l 



n 



+ 



27ri z — \ 



\Q{z) 



Case a = 1 



Let 



P(z) 



-z^[l + {z-l)E] 



ihz - 1) 
Then if Zk is the root such that Zfc — ^ 1 as C — i> with for small C 



zu = l+ Ce'^^'Pizk 



and if h{z) is analytic near z = 1, one has 



/i(^fc) = f Y~^^^ 



1- 



C'e2-^P(zfc)l " P'(z)/[A'/P(2)] 



;z - 1 - [Ce2^-'=P(zfe) 



Then summing over the roots < k < M — 1 leads to 
^"^ °° 1 r dz 



k=0 n=\ "' •'^ ^''^ 



(^)- 



[P(z)]' 



(A16) 



(A17) 
(A18) 



(A19) 



(A20) 



(A21) 



(A22) 



(A23) 



{z - 1)"*^ 
Therefore the equations for case a = 1 are given by exactly the same expressions as the case 



a 7^ 1 (^,|J,^) with the replacement 



+ 
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Appendix B: Asymptotic Evaluation of Velocity and Diffusion Constant 

Evaluation of the Velocity 

For 1/a < Zc and 1/6 > Zc 

In this phase (^) is dominated by Slm the saddle point contribution 

XL,M = ^L,A/ + 0(Sz.,M/i:) (Bl) 

where 

^'^^ " v/2^(6 + p-l)(« + p-l)(^e-l)^- ^ ^ 

However the leading contributions to Xi^^m cancel in (|y). To evaluate the next leading 
contribution we write Zi^^m as a double integral using (|38D 

7 = I ^^ ^ I ^^ ^ ^^ z^ ^^~^^^^~4 p.->.1 

^'*^ / 27ii {bz - l){az -1)1 2m {h~z - l){a~z - I) {z - 1)^^ {~z - ly^ 1 ^ ' 

The double integral is dominated by the saddle point z = z = Zc = 1/(1 — p) and can be 
evaluated to be 

(B4) 







pSIm 
(l-p)L 


(|i2|) one obtains v = 


= l/Zc - 


- {zc — '^)/ Zc SO that 
f = 1 - 2p 



(B5) 
as in (^. 
For l/a < Zc and 1/6 < Zc 



As explained in section 6, in this phase the dominant contribution to (]38|) is the pole at 1/6 
and the saddle point is sub dominant: 

Xl,m = Bl,m + Sl,m + 0{Sl,m/L) (B6) 

where B represents the contribution of the clockwise contour around the pole z = 1/b 

^^•^ = 6- a (1-6)^^6^-*^^- ^^^^ 

Therefore, from (^) 

(1 - bzcf 

Zl,m — 7 BlmSl,m (B8) 

bzc 
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and from (^21) 

v = l- (3- p (B9) 

as in (^. 

For 1/a > Zc and 1/6 > Zc 

In this phase the pole at 1/a is the dominant contribution to (|38D and the saddle point is 

sub dominant 

Xl,m = Alm + Si^^M + 0(5i,M/^) (BIO) 

where v42,,Af is given by 

1 1 

Zl^m may be determined by symmetry considerations from the previous phase: under in- 
terchange of a and b ( PB| ) becomes ( [BIOD and (BH) becomes 



and we find 

as in (p!0|). 

For 1/a > Zc and 1/6 < Zc 
In this phase 

so that 

and 

as in (PH). 



'^L,M — ^^AlmSl.m (B12) 

azr. 



V = a — p (B13) 



Xl,m ^ Al,m + Bl,m + 0{Sl,m) (B14) 

Zl,m ^ ^ ^ ^ A^,M^L,M (B15) 

ba 

v = a-P (B16) 



Evaluation of Diffusion Constant 



In order to compute the diffusion constant given by (^ 



^ = ^r^f^L,M (B17) 



'LM 



where 

UrM 



W2L,2M-iZl-IM-1 + W^2L- 



l,2M-2\^L-l,M 



JL-1,M-1, 



W2L~2,2M-2Zl,M ■ (B18) 



we need first to evaluate the asymptotics of Ul,m- Using the integral definitions (|4^JB3|) we 
may write ( [B18|) as 



U, 



L,M 



dw 



w 



2L 



[1 + E{0){w - 1)]^ r dz 



1 



X 



2m {w - If^ {bw - l)^{aw - l)^ / 27ri {z - 1)^ {bz - l){az - 1) 
dz z^ {z — \){z — z) w — lw — zw — z 



P 



(B19) 



2m {z - ly^bz - l){az - 1) w^ z 

For 1/a > Zc and 1/6 > Zc 

In this phase all integrations are dominated by their saddle points. However care is required 
to correctly identify the first non- vanishing contribution to ( pi9|) . This is done most sys- 
tematically by considering the scaling of the large deviation function and we carry this out 

in section 7 where we show 

(Lvrp(l-p))V2 



A 



(B20) 



For 1/a < Zc and 1/6 < Zc 



To evaluate ( |B19| ) we carry out the integrals in sequence. In the first integral over z we 
keep an apparently subdominant term (proportional to Sl,m) as well the term proportional 
to Bl,m, because when we integrate over z we find both terms give leading contributions 
proportional to Sl,mBl,m i-e. the dominant contributions to the triple integral come from 
w at the pole, one of z, z at the pole and the other at the saddlepoint: 



U. 



L,M 



dw 



w 



2L-2 



[l + E{0){w-l)f r dz 



1 



2m {w - l)2M-i (py^ _ i)2(q,^ _ 1)2 J 2m {z - 1)^ {bz - l){az - 1) 

{b-l){bw-l){w- z){bz-l) {zc-l){w- Zc){w- z){zc- z) 

^L,M — r Ol,M' 



bz 



z'^z 



dw w^^-^ [l + E{Q){w-l)f 
2m {w - l)2A^-i {bw - lY{aw - If^^^^' ^'^ 
{b - l){bw - l){w ~ zc){bzc-l) {zc-l){w - Zc){bw - l){bzc - 1) 



Sr,krB 



L,A/-DL,A/- 



Slm^l.m 



{bzc 



bz, 

I? 



bzl 



bzl 



dw 



w 



2L-2 



[l + ^(0)(w;-l)]2 {w-Zc 



2Tii {w — 1) 



2A/-1 



'l-b){bZc 



{bw — l){aw 

r \ n 2 



1 + ^(0) 



(B21) 



In this phase the behaviour (B6) and the form of (B7,B^) imply that -E'(O) given by (p7|) 
becomes 



E{Q) 



1 + 



:i - bZr) S, 



[l-b) Bl 



L,M 



M 



(B22) 
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Therefore 



and (|B17| ) along with 



Ul,m 
6D and h = 1 

A 



St AfBr 



{hz, - If 



/3, yields 
p-(3 



(B23) 



(B24) 



as m i^. 

For l/a > Zc and 1/6 > Zc 

For this phase the evaluation of the integral (|B19| ) is very similar to that outlined above 
for the previous phase, with a replacing h. In carrying final integral over w a factor of —1 
is introduced due to the opposite directions of the integral around l/a and 1/6. Therefore 
one obtains the diffusion constant by interchanging (3 and 1 — a in (B24) and multiplying 
by -1 

A ^ ^il^^ . (B25) 



i — p — a 
as in (|T0|). 

For l/a > Zc and 1/6 < Zc 

In this phase it turns out that the dominant contributions to Ul,m come from one of z, z at 
the pole 1/6 and the other at the pole 1/a and w at either of the two poles. Carrying out 
the integrals in ( [B19| ) in sequence gives 



U 



L,M 



dw 



w 



2L-2 



[1 + E{0){w - 1)]'^ r dz 



\M 



2m {w - l)2M-i (5^ _ i)2(^aw - If J 2m {z - 1 

(a — l){w — z){aw — 1) {b — l){w — z){bw — 1) 

^L,M n ~ 1" Bl,M 



-^l.mBlm- 



a{bz — l)z 
a)"^ r dw 



ab 



b{az — l)z 

y,2L-2 [l + ^(0)(w;-l)]2 



2m (w 



Al,mBlm 



ib-afil-b) 



a 



l)2A/-i (bw~l){aw-l) 
\ n2 



1 + ^(0) 



1 



~^'^L,M-^^^'^' 



{b-af{l-a) 



m) 



1 — a 



a 



-,2 



(B26) 



First consider B^^m ^ Al,m- Then due to the form of X^^m ( 1B14| ) in this phase one has 



E{0) 



1 + 



(a - 6) Al,m 



Both terms in 



U, 



1-6 

contribute and one obtains 
{b-af 



1 - b)a Bl, 



M 



L,M 



AI,,Bl,m^,^^^_^^, 



[6(l-6)+a(l-a)] 



(B27) 



(B28) 
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Then (pT^ ), (|BT5| ) and (|BT7|) imply 



A ^ ^'^-^'+°'^-°' (B29) 

1 ~ P — a 



as in (pup. 

In the case where Al^m ^ Bl^m 

E{0) - 



a 



1 — a 



a) Bl. 



M 



6(1 - a) Al, 



M 



(B30) 



however it turns out that one obtains the same expression for the diffusion constant (P29|) . 
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